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CN ■ We use Fuchsian methods to show that, for any two dimensional manifold E 2 , there is a 

' large family of U(l) symmetric solutions of the vacuum Einstein equations on the manifold 

, S x S 1 x K, each of which has AVTD behavior in the neighborhood of its singularity. 

; 1 Introduction 

7— I . 

■ In l]. Isenberg and Moncrief obtain a family of U(l) symmetric cosmological solutions of the 

vacuum Einstein equations which have asymptotic velocity term dominated ( "AVTD" ) behavior 
in the neighborhood of their singularities. Such behavior, which has been studied in a number of 
classes of spacetimes (See, for example 0, 0! and 0.)> has been found to be useful for the study 
of the strong cosmic censorship conjecture and similar questions [3]. Note, however, that little 
' is known about the long time evolution of these spacetimes away from the singularity; indeed, 

the proof PP that there are U(l) symmetric solutions with AVTD behavior relies on the use 
| of Fuchsian methods, which only guarantee existence of the solutions in a neighborhood of the 

O -1 singular region. 

In H 0, Choquet-Bruhat and Moncrief obtain a family of U(l) symmetric cosmological 
solutions of the vacuum Einstein equations which expand for an infinite amount of proper time, 
and which asymptotically approach (in the expanding direction) a flat spacetime. In this case, 
little is known about the behavior of these spacetimes in the contracting direction; a singularity 
generally forms, but it is not known whether the gravitational fields show AVTD behavior in the 
neighborhood of the singularity. 

Are there solutions which combine these two behaviors: AVTD behavior near the singularity 
in the contracting direction, and global existence with asymptotic approach to a flat spacetime in 
the expanding direction? If one attempts to show that such solutions exist by somehow combining 
the results of [Hj with those of p^ cited above, one immediately runs into a problem: The AVTD 
results of pQ are formulated for spacetimes that are topologically T 3 xR, while the global existence 
and stability results of [5] work with spacetimes that are topologically S 1 x S 2 x R, where E 2 is a 
surface of genus greater than one. The question then becomes whether either of these results can 
be adapted to alternative spacetime topologies. In this work we show that for the AVTD results 
this is indeed the case. 

We show, using a global in space formulation which is independent of the topology, that for 
any surface S 2 , there is a family of U(l) symmetric vacuum spacetimes on S 1 x E 2 x K which 
have AVTD behavior in a neighborhood of their singularities. To obtain this result, we have had 
to add a restriction beyond the polarization (or half-polarization) condition which is imposed in 
[P. This restriction, which we spell out in section 4, controls the asymptotic (in time) conformal 
class of the spatial geometry transverse to the orbits of the U(l) isometry group. As a result, one 
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can study (globally in space) the dynamics of the geometry near the singularity, regardless of the 
topology of E 2 . 

The proof of our results here, like that used to prove the main results in , relies on Fuchsian 
methods. We carry out the proof in sections 4 and 5, after describing the spacetimes of interest in 
section 2, and discussing the velocity term dominated equations and solutions for these spacetimes 
in section 3. We state our main theorem in section 5, and then make concluding remarks in section 
6 regarding the work that needs to be done if we hope to find solutions with both AVTD behavior 
in the contracting direction and global existence in the expanding direction. 



2 Polarized U(l) Symmetric Vacuum Spacetimes 

It has long been known that the vacuum Einstein equation system on a spacctime ( 4 TVf, 4 g) with 
a U(l) isometry group acting spatially is equivalent to a coupled system of the Einstein equations 
and a wave map on the reduced spacetime ( 3 M, 3 g). Further, if we assume that the spacctime 
( 4 M, 4 g) is polarized in the sense that the 3-planes orthogonal to the orbits of the isometry 
group C7(l) are surface-forming, then it follows that the wave map reduces to a (real) scalar field 
satisfying the wave equation on ( 3 M, 3 g). 

To be more explicit, we choose a smooth surface E 2 (of arbitrary genus), we set M = 
S 1 x E 2 x R, we choose 8 as the coordinate for the circle (the orbits of the U(l) action), and we 
write the metric 4 g in the form 

* g = e- 2 *( 3 g) + e 2 *d6 2 , (1) 

where 3 g is a Lorentz signature metric on E 2 x K, and <fi is a scalar field on the same (three 
dimensional) spacetime manifold (both independent of 9). The vacuum Einstein equations on 
4 M now reduce to the system 

R a p[ 3 g] = 2d a <j>dp4> (2) 

D-D(f> = 0, (3) 

where D is the (three dimensional) spacetime covariant derivative corresponding to 3 g, D ■ D 
is the corresponding D'Alembertian (wave) operator, and the Greek indices a,f3, etc, take the 
values {0,1,2}. 

We wish to reformulate the Einstein system (HJ-© in 2 + 1 dynamical form. We first rewrite 
the 2 + 1 spacetime metric in the form 

3 g = -N 2 dt 2 + 7 , (4) 

where TV denotes the lapse, 7 is a (time-dependent) Riemannian metric on E 2 , and we have 
(without loss of generality) assumed that the shift is zero. The dynamical variables now consist 
of the metric 7 together with the extrinsic curvature 

Kb = -j^dtlab (5) 

(Here and below the Latin indices a,b, etc, take the values {1,2}; they denote the indices for a 
coordinate basis on the 2-surfaces E 2 .), and the scalar field (f> together with 

P = JjM- (6) 

The lapse is not a dynamical variable; its evolution may be freely chosen, and corresponds to a 
gauge or coordinate choice. 
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It is useful for later purposes to replace the variables (-/V, 7) by the equivalent set (A,er), with 
the transformation between these two sets of variables given by 



N = e x (7) 

and 

lab = e X (T a b- (8) 

We note that this replacement of (N, 7) by (A, a) does not imply any restrictions on the fields, 
and does not entail any choice of gauge. We also note that in terms of these new variables, we 
have 

Kb = -^{d t <Jab + <7abd t X), (9) 

t = l ab k ab = -e-\d t X + (10) 
where for convenience we use the definition ip '■= u ab dt<J a b, and 

p = e- x d t cj ) . (11) 

As usual for the n + 1 dynamical formulation of Einstein's theory (here, n = 2), the field 
equations split into constraints and evolution equations. The three constraint equations for this 
reduced system take the form 

R(l) - k ab k ab + t 2 = 2p 2 + 2 7 ab 3 a </)d 6 </> (12) 

and 

V b k b - d a r = -2 P d a & (13) 
while the evolution equations can be written as 

d t k ab = NR ab - 2Nk ac k c b + Nrk ab - V a d b N - 2N8 a <j)d b ^ (14) 
dtp = -\p4> ~ pd t X + e V(V a V^ + V b <RA), (15) 



together with 
and 



daab = -2e A k ab (16) 



dt<f> = e x p. (17) 

Here V denotes the covariant derivative compatible with the metric 7, R ab is the corresponding 
Ricci curvature, and R := j ab R ab is the corresponding scalar curvature. 

We have not yet chosen a gauge condition; we do that below. We note, however g^T^g = 
7}'4' e ~ 2X > so the system is in time harmonic gauge if and only if tp = 0. 

As noted in the Introduction, in this paper we impose an extra restriction on the class of 
U(l) symmetric spacetimes we are examining: We require asymptotic (in time) control of the 
conformal class of the two-dimensional spatial geometry transverse to the U(l) group orbits. 
In section 0] below, we make this requirement (which we call the "asymptotic conformal class 
condition") mathematically precise. 



3 



3 VTD Equations and Solutions 



A spacetime has AVTD behavior if, as one approaches the singularity (chosen here to be at 
t — > oo), the fields approach a solution of the VTD equations. One obtains the VTD equations 
from the Einstein equations by dropping all terms involving space derivatives, except in the 

momentum constraint (|13[) . Thus, if we use the notation ( A / a bi'j ) ,k a b,'p) < ► {&ab, A, 4>, k ab ,p) to 

denote these fields appearing in the VTD equations, we have the constraints 

- k ab k ab + t 2 = 2f (18) 

and 

V b k b a - d a r = -2pd a $, (19) 

and we have the evolution equations 

dtlab = -2e x k ab (20) 

and 

d t k ab = -2Nk ac k c b + Nfk ab , (21) 

together with 

d t 4> = e' x p. (22) 

and 

dtp = -jP4> -pd t \- (23) 
Note that to solve the VTD equations, we find it useful to rewrite <|21ll in mixed index form: 

d t k b a = Nrk b a . (24) 

To enforce the asymptotic conformal class condition on these VTD solutions, we set the 
quantity a ab equal to a time-independent field a ab on the 2 dimensional manifold E 2 : 

Oab = (Jab- (25) 

(We shall generally use the tilde to denote time-independent quantities.) It follows immediately 
from its definition that "0 = 0, and it follows from (JSJ and from © that all of the time dependence 
in and k ab is contained in A. Consequently, from (|24(l . we deduce 

8f t X = 0, (26) 

which has the general solution 

X = X — vt, (27) 

with A and v a pair of arbitrary time- independent functions on S. The general solutions for j ab 
and k ab now follow. We have, for example, k £ = ^e~ x S^v 

It remains to solve for </> and p. With t/j = 0, we may combine (|17|l and Ij23(l to derive 

dU = 0, (28) 

from which we obtain the general solution 

4> = 4>-wt (29) 
for the time independent functions <f> and w. It follows that p = —e~ x w. 
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What about the VTD constraint equations? If the solutions for k ab and f and p are substituted 
into the VTD Hamiltonian constraint l|18|) . then it reads 

e- 2 \ l -tf-2w*), (30) 

which implies that | iZ> | = \\v\. This reduces, up to signs of v and w the set of free functions in 
the VTD solution to v, </>, A, and a ab . 

A similar substitution of the VTD solution quantities into the momentum constraint H13(l 
produces the equation 

- v(d b X - d b vt) + d b v + Awd b (4) - wt) = (31) 

(Here we need the 2 + 1 identity 3 r^( 7 ) = T c ab (a) + |(^9 Q A + S c a d b X - a cd a ab d d X).) We treat the 
case w = -xV, then the momentum constraint can be written in the form 

d b [ve { - U2 ^} = 0, (32) 

which permits us to solve for v, and thereby eliminate one more of the free functions in the 
asymptotic data; we are left with just 0, A and a ab . 



4 AVTD Expansion and the Fuchsian Analysis 

The procedure for using Fuchsian PDE methods to prove the existence of various families of 
solutions of Einstein's equations characterized by AVTD behavior near a spacclikc singularity is 
by now fairly familiar [2] PP 0j . One starts by writing out expressions for each of the fields as the 
sum of a solution of the VTD equations plus a remainder function multiplied by an exponential 
decay factor. This is done with no loss in generality, since the remainder fields are completely 
general. One then substitutes these expressions into the Einstein equations (first order form), 
and verifies that (perhaps with certain restrictions assumed on the free functions in the VTD 
solution) the remainder fields satisfy a PDE system in Fuchsian form [5] . It follows that for each 
of the allowed solutions of the VTD equations, there is a decaying solution for the remainder fields 
in the neighborhood of the singularity; consequently there is a solution of the Einstein equations 
which asymptotically approaches the given solution of the VTD equations. 

We find that this analysis is a bit simpler if we make a slight change in the choice of the field 
variables we use to describe our U(l) symmetric spacetimes: we replace a ab by its inverse a ab , 
and we replace k ab by k b (index raised by the metric inverse r y ab ). Then we write the expansion 
formulas for the fields as follows: 

a ab = a ab + e-'° t 8a ab , (33) 

A = A - vt + e- e ^SX, (34) 
cf) = 4> - wt + e- e ^S4>, (35) 

k b a = e- x (^v6 b a + e-^Sk b a ), (36) 

p = e- x (-fv + e- e " t dp). (37) 

Here e a , e\, e^, e&, and e p are all positive numbers to be chosen later, while 5a ab , SX, S(j>, Sk b , and 
dp are the remainder fields. 

At this stage, it is useful to make a choice of gauge. Recalling that a choice of the lapse and 
shift largely control the gauge, and recalling Q, we set SX = 0. That is, for a specified choice of 
the VTD data, we set 

A = A = A - vt. (38) 
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For later purposes, we note that this gauge condition i|38|) is equivalent to the condition 

±i/> + e-^SkZ = 0. (39) 

We recall that the shift has been set to zero. 

In checking whether the Einstein equations together with the above expressions for the initial 
data fields result in a Fuchsian PDE system for the remainder fields, we recall that to do this we 
must cast the system into first order form. Thus we need to add extra field variables, representing 
the spatial derivatives of those fields-0 and a ab - whose second derivatives appear in the Einstein 
equations H12fl - I|15|) . To take care of the first of these, we introduce the new variable ^>r c i, and 
a) everywhere first or second derivatives of 4> appear in the Einstein evolution or constraint 
equations, we express them in terms of 0r c i and its first derivatives; b) we assign to </>r c i the 
evolution equation 

dt<t> [c ] = d c (e x p); (40) 
and c) we assign to 4>[ c ] the expansion formula 

[c] =0 c 0-0 c iW + e-V*ty [c] . (41) 

Note that the expansion formula l|41|l guarantees that if the remainder terms are bounded for 
both d c 4> and <j)M, then d c 4> and 4>[c] agree asymptotically. 

Similarly, we also introduce the new variable cr^j, and a) everywhere first or second derivatives 
of a ab appear in the Einstein evolution or constraint equations, we express them in terms of 
tTjJ-j and covariant derivatives V m of of| (Here V m denotes the Lcvi-Civita covariant derivative 
corresponding to the asymptotic metric a a b)', b) we assign to er^j the evolution equation 

dtoft = V c (d t a ab ) = (-d c v)a ab va$ + 2e x [k h m a^ + a am V c k b m + (d c \)a am k h m ], (42) 

where we calculate dtd ab using a ab = e x -f ab , dtjab = — 2e A fc a f,, and (|38|l : and c) we assign to cr?j 
the expansion formula 

^ = e-^Safl (43) 

Again, these formulas are chosen to result in asymptotic agreement of cr^ and V c er af ', presuming 
that the remainder terms are bounded. 

We are now ready to verify that the induced equations for the vector quantity 

U := (Sa ab , d>, Ski Sp, 6v$, 6<j> [c] ) (44) 

take the (first order) Fuchsian form 

d t U -LU = e~ flt F(x,t,U,d a U), (45) 

where L is a linear operator which is independent of t and has negative eigenvalues, fi is a 
positive number, and F is a function which is uniformly Lipschitzian in U and V?7 so long as 
those quantities are uniformly bounded, and is analytic in x and t. We do this verification for 
each of the components of U separately, since as we shall see, the operator L is diagonal. 

4.1 Equations for 5a ab and (Sof-j 

To derive the evolution equation for 6a ab we work with the equation dtj ab = 2e x k ab . For the 
right hand side, we use I® and J33J| to calculate 2e x k ab = e" A (ucr ah + 2e~ 2tkt a ac 8k b ) . For 
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the left hand side, we use ggj and (J33J) to derive <9 t 7 Qb = <9 t (e~ A cr Qb ) = e - x (i a ab + d t a ab ) = 
e~ x (vcr ab — e <T e~ <L " t 5a ab + e~ € ^ t d t Sa ab ). Then combining, we obtain 

d t Sa ab - e a Sa ab = 2e (t °-^ a am 6k b m (46) 

This equation matches the requisite Fuchsian form so long as one chooses > e CT > 0. 

To obtain the evolution equation for 5<r^ we substitute the various expansion formulas into 
(I42|l . We obtain, after some work, 

d t S<T$ - e a ,G\\ = 2e^'-^ t V c {a bm ~5k a m ). (47) 

We have Fuchsian form here so long as one chooses e k > e a > > 0. 

4.2 Equation for 5k b a 

Our derivation of the form for the evolution equation for Sk b is based on the evolution equation 
for k b , which is as follows 

d t k b = NR b + Nrk b - N\7 b \7 a N - 2N 1 bm ( j >[m] cj >[a] . (48) 

Using the expansion formulas, we now look at each of the terms in l|48(l : We first calculate the 
time derivative term on the left hand side. From (|3fc>|) . we obtain 

dtk h a = e~ x [\v 2 5 b a + (v- e k )e-^Sk b a + e -^d t 6k b a }. (49) 
We next calculate the Nrk term 

Nrk b a = er x [\f5 b a + i e -^Sk b a + V^^fc™^ + e^5k^8k b a ] . (50) 

Note that in this expression, as in (05J, all but the first term has a factor e~ £fct (The first terms 
in these two expressions cancel). In the expressions we calculate below, no such factor appears. 

For the hessian of N, we find 

V b d a N = e- x a bm S7 m d a e x = a bm [d m Xd a X + 8 m d a \ - T d am ( 7 )d d X\, (51) 

where rf m (7) are the Christoffel symbols calculated from the metric 7& c . We could now expand 
this expression using the expansion formulas for a ab , a? b y and A. However, the resulting formula is 

cumbersome, and all we really need to know is that it is analytic in 5a and ^c^j, and that it has 
coefficients with at most polynomial growth in time. A similar analysis holds for the remaining 
two terms. For NR^, since we are working on a two dimensional manifold, we have 

NR b a = e x R b a = \e x 5 b a R(l) = ^ b a [R(c7) - A a X], (52) 

where A a is the laplacian corresponding to <r, and R{<j) is its scalar curvature. This expression, 
if expanded out, is again analytic in 8a ab and Sa^, and V d^ ^ a s well, with coefficients growing 
no faster than polynomial in time. Finally for — 2N-f bm (j)i rn \^>i a \, we readily determine that this is 
analytic in 8a mb and (S0[ c j , with the same sort of coefficients. 

Putting all of these expressions together, we have 
d t 5k b a - e k Sk b - \v5 h a 8k™ = e^A^fe* (53) 
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where F b is a function of the indicated variables, and is polynomial in 5a ab , Sa^, ddScr?^!, S(f>, and 
<5</>[ c ] with bounded coefficients depending on x and t. Note that in deriving (|53[) . we have used 
the expansion e (A+efc * ) = e (efc-D) *e\ 

To verify that this equation is of Fuchsian form, we find it useful to split it into its trace and 
trace-free pieces. For the trace piece, we have 

d t 5kl - (e k + i)5k b = e-^(SkZ) 2 + e^-^F b (x, t, Sa a \ Saf^ d d 8o$, 5<j>, S<f> [c] ). (54) 

Clearly this is of Fuchsian form so long as v(x) > 6k > 0. A similar observation holds for the 
trace-free piece. 

4.3 Equations for 5p, 5<p and 5(j)^ 

We derive the evolution equation for the "scalar field" variable dp (recall that </>, (j>M , and p 
are really geometric variables) from the field equation (|15fl (with appropriate replacement of 
derivatives of 4> and a ab by </>[ c ] and cr?X) . To carry this out, we calculate each term in (|15|l using 
the expansion formulas for the relevant field variables, along with the imposed gauge condition 
(EH- 

We first calculate the term on the left hand side of equation (|15|l , dtp: 

dtp = dt[e- x (-w + er^dp)] = e - x [v(-w + e^Sp) + e- ept (-e p Sp + d t Sp)]. (55) 
Next we obtain 

- \p<P -pd t X = e- x [-w + e- e »*6p][v + e^SkZ}- (56) 

and 

eV b V a [b] - cr ab (d a <f> [b] - 3 r^ [c] ) = a ab d a <j> [b] - a ab T c ab {a)4> [c] , (57) 

where 3 r^ b denotes the 2 + 1 dimensional spacetime Christoffel symbol, and T ab denotes its 2 
dimensional counterpart. These are related by the expression 3 T ab = T c ab + ^(S b d a X + S^dbX — 
a cd a a bddX); however we note that in the calculation done in (|57|) . the A terms cancel. As with 
the hessian and curvature terms in the dt5k b a calculation, the full expansion for e x j ab V a <P[b] is 
rather cumbersome. Again, it suffices to note that when we carry out this expansion, we obtain 
an expression which is analytic in 5a ab , Sa^, 8(j), 54>[ c ], and ddS(j)[ c ], with coefficients that have at 
most polynomial growth in time. 

Combining all of these terms, we obtain 

d t Sp - e p 5p = -e^-^tiSkZ + e~ ekt S P SkZ + e (e »- D) *F(x, t, 5a ab , Sa$, 6(f>, 8^ [c] ,d d S<j) [ch 5p), 

(58) 

where F is a function of the type just described. We see that (|58|l is in Fuchsian form so long as 
£fe > £p > 0, and v > e p . 

The remaining evolution equations are very easily checked. From the equation dt4> = e x p, we 
obtain 

d t 5(j) - e^Scj) = e^f-^Sp, (59) 
which is clearly of Fuchsian form so long as e p > > 0. Then from dt4>[ c ] = 9 c (e x p) we derive 

cW [c] - e^ [c] = e^t'-^dcSp, (60) 
which is Fuchsian so long as e p > e^/ > 0. 
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4.4 Applying the Fuchsian Theorem 

Taken together, the results of sections 4.1-4.3 show that, for any choice of the positive constants 
£ := (e a , efc, e a i, e<p, e p , e^i) and for any choice of the asymptotic data A := (cr ab , <fi, A, v, w) satis- 
fying the conditions 

v(x) > € k > Max{e a ,e a >,e p } (61) 

and 

e p > Max{e^ e^} (62) 

the PDE system consisting of lj4U)> . (|4T|) . lEHJ; ESt ; an d (09 satisfied by the components of 
U = (5<r ab , Stfi, Sk b , Sp, 5a^Y S(f>[a]) is of Fuchsian form. The following result is then a consequence 
of standard theorems regarding Fuchsian systems [S] 

Proposition 1 For any choice of the positive constants £ and analytical asymptotic data A 
which satisfy the conditions i)67)) - there is a unique analytic solution U{t) of the equations 
\4b\ , {4'$ , \5tl\) . }58\) . }5y\) . and \6U\) for sufficiently large t. Further, this solution has the limit 
lim^oo U(t) = 0. 

5 AVTD Solutions of the Einstein Equations 

Since we have derived the Fuchsian system (QUI, igJJ, (JSHJ, (JSSJl, $5$, (jHOJ) by substituting the 
expansion formulas for (er ab , erg, k£, A, <fi, 4>^ c ^p) into the Einstein evolution equations, we should 
be able to infer from Proposition 1 conclusions regarding solutions of Einstein's equations. To 
do this, we first need to verify that the unique solution U whose existence is guaranteed by 
Proposition 1 satisfies d c (f> = 4>[ c ] an d V c cr oh = 

To make this verification, we use the evolution equations (|42|l and l|40(l to show that dt(4>[ c ] — 
d c 4>) = and dt(o-^ — V c er afc ) = 0. Hence the quantities </>[ c ] — d c <fi and cr^ — \7 c a ab are independent 
of time. Next we verify that lim t ^ 00 ((/)[ c ] —d c 4>) — lim t _ >00 ([9 c (/) — d c wt J r e~' L <f>' t 54>ic]\ — [d c (f)~d c wt+ 
e~ e 'i' t d c 6<p}) = lim t ^ 00 (e~ e i>' t S(f>- e~ e **<9 c <5</>]) = 0. Similarly we calculate lim t ^ 00 (a-?| - V c ct q6 ) = 

lim^o^e^'Mcrg - e- e "*V c Sa ab ) = 0. It then follows that d c <p(t) = (f> [c] (t) and V c a ab (t) = erg 
for all time t during which the solutions exist, and so we have proven that we have a family- 
parametrized by the asymptotic data A -of solutions (jabit), fc c d, 4>(t),p(t), N(t)) of the Einstein 
evolution equations for sufficiently large t, and further, the solutions are AVTD in a neighborhood 
of the singularity at t — > oo. 

It remains to show that, for certain choices of the asymptotic data A, these solutions of the 
Einstein evolution equations also satisfy the Einstein constraints, and hence the full Einstein 
system. We argue this via the following chain of observations and lemmas. 

We use the label S_A(t) to denote the unique (AVTD) solution of the Einstein evolution 
equations which has asymptotic data A, and we use the notation C M :— (Co,C a ) = (— h[R — 
k b k a + t 2 - 2p 2 - 2|V0| 2 ], -[Vbfck - V a r + 2pV a 0]) to denote the three constraint functions. 
We first note the following claim (which one proves by substituting the expansion forms of the 
dynamical variables into the constraint quantities C M ): 

Lemma 1 For any choice of A, the corresponding solution 1S4 (t) satisfies the limit equation 

lim C^(S A (t)) = C^A), (63) 

where denotes the VTD constraints \3U\) - \31)) . 
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For a general choice of A, C^A) is nonzero. However, as shown in section 3, the VTD 
constraints are satisfied if and only if w = iv, and v = ce^ x ~ 2 ^' for any constant c. We denote 
by Ao those sets of data for which this is true. 

We now seek a first order system of evolution equations for the constraint functions C M . 
While one could derive these using the chain rule and the evolution equations for the dynamical 
variables, for present purposes we obtain the system in its most useful form by relying on the 
contracted Bianchi identity, which takes the form D V G V ^ = 0, where G ufl is the 2 + 1 spacetime 
Einstein tensor. Combining this with the wave equation D ■ Dtp = for 0, we have 

D v % = 0, (64) 

where := G% — T" so 0^ = corresponds to the Einstein equations. In particular, we have 
Oq = Co, and 0° = e~ A C a where Co and C a are the three constraint functions, written explicitly 
above. Expanding out (|64|l . and using the 2 + 1 decomposition expressions for the connection 
coefficients, together with some of evolution equations for the dynamical variables, we obtain 

d t C - 2e A rC - 7 ah V a (e A C 6 ) + e V fc (V a A)C b (65) 

and 

8 t C a - e\C a = e A V Q C + 2e A (V a A)C . (66) 
Setting Eq := e 2A Co and E a := e A C a , we may rewrite these equations in the form 

d t E - 2vE - 2e A r^ = e x a ab \7 a E b + <7 o6 e A (V a A).E a (67) 

and 

d t E a - iE a - e\E a = V a E . (68) 

We immediately note that, for any choice (solutions or not) of the functions A, a ab ,T and 
v which appear in the coefficients of H67f) - i|68|) . these equations are linear and homogeneous in 
(E ,E a ). Hence (E , E a ) = (0, 0) is always a solution. We want to show that for any data Aq for 
which the asymptotic condition C M (^4o) = holds, C^(Sj^ ) — necessarily holds. The key step 
towards that goal is to verify 

Lemma 2 // the coefficients for the system j6'7| )- f6'<^)) are constructed from fields Sa for which 
the asymptotic condition C^{Aq) = holds, then the system j6'7| )- f6'<S)) is Fuchsian 

Proof Since the coefficients of these equations are constructed from solutions of the evolution 
equations, we may use the expansions and estimates derived in section 4. Hence, we have e A r = 
v+e~ ekt 5k™, from which we determine that the left hand side of (|FT7j> is dtEo—4:vEo — 2e~ €kt Sk™Eo 
and the left hand side of is dtE a — 2vE a — e~ Skt 5k™E a . Removing the last term of each of 
these expressions (instead placing each on the right hand side of the respective equations), and 
recalling the inequality (|61|l . we see that the remaining expressions dfEo — AvEq and dfE a — 2vE a 
are clearly of proper Fuchsian form. Recalling that e A = e x ~ vt , and noting the boundedness 
of a ab and of e~ £kt 5k™, we see that the right hand side of (|67|l together with the added term 
4e~ ekt 8k™Eo is of proper form as well. It remains to check the right hand side of l|68|) . together 
with its added term 2e~ tkt b~k™E a . 

The added term is again clearly no trouble. As for the rest, writing out E in terms of 
dynamical variables, we have 

E Q = e 2A C = -\e 2X [Rd) - k a b k b a + r 2 - 2p 2 } (69) 
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It follows from Proposition 1 that the quantity in the bracket in (|69[) is bounded for large t. Hence 
we can write 

\W a E \ < e-^F(x,t), (70) 

for some positive constant /x and some bounded function F(x,t). This establishes that the right 
hand side, and therefore all of, (|68|) is of Fuchsian form, so we have proven the Lemma. 

Since the evolution equations for the constraint functions are (for appropriate asymptotic 
data) of Fuchsian form, it follows from the standard theorems that for such data, there is a 
unique solution which decays asymptotically. We already know, however, that (Co,C a ) = (0,0) 
is a solution. Therefore, (Co,C a ) = (0,0) is the only solution. We conclude that the constraints 
are satisfied by the fields Sa , and we summarize our results as follows: 

Theorem 1 Let E 2 be any two dimensional manifold. For any set of asymptotic data Aq on E 2 
which satisfies the asymptotic constraints and the asymptotic conformal class condition, there is a 
unique polarized U(l) symmetric solution of the vacuum Einstein equations on S 1 x E x R which 
has AVTD behavior in the neighborhood of its singularity. 

6 Conclusions 

The existence of a large family of 1/(1) symmetric solutions with AVTD behavior has been es- 
tablished in earlier work pQ. The significance of this work is that we now know that for any two 
dimensional manifold E, there is such a family of AVTD solutions on E x S 1 x R. 

As noted in the Introduction, the goal of this work is to be able to produce solutions with 
AVTD behavior in a neighborhood of the singularity, and with global (proper time) existence 
evolving away from the singularity. One possible way to do this is to show that there is an AVTD 
solution of the sort discussed here, which admits a constant mean curvature (in the two dimen- 
sional sense) surface for which an appropriate energy for the data on these surfaces satisfies a 
smallness condition. It is not yet clear whether or not this can be done. 
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